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Abstract 


In this work, the characteristic basis function method is adapted to analyze the problem of scattering 
from multiple and multi-scale impenetrable targets in conjunction with the discontinuous Galerkin 
method, and the monopolar RWG functions are chosen as the basis functions. This method enables 
us to analyze multiple and multi-scale targets using nonconforming discretizations. The use of the 
CBFs helps reduce the size of the impedance matrix of the associated method of moment significantly, 
enabling us to employ direct solvers as opposed to iterative solvers. The process of generating the 
reduced matrix is naturally parallel, and the reduced matrix is well conditioned, which obviates the 
need for pre-conditioning. In addition, the adaptive cross approximation algorithm is implemented to 
reduce the complexity of the computation. Numerical results are included to demonstrate the 
accuracy and efficiency of the present approach when analyzing scattering from multiple and multi- 
scale targets using nonconforming discretizations. 


Index Terms: Characteristic Basis Function Method (CBFM), Integral Equation Discontinuous Galerkin 
Method (IEDG), Adaptive Cross Approximation (ACA), Method of Moments (MOM), Size Reduction, 
Multiple Scale. 
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What's this? 


«¢ The analysis of electromagnetic scattering 
from targets dividing into parts,is to 
numerically analyze the arbitrary-shaped 
target whose component parts divide from 
each otherin practical engineering, like 
rockets and missiles. 

» 





Why do this? 


Reduce the unknowns and 
memory consumption, accelerate 
the computation. 


Use only one mesh to analyze the 
target dividing into parts. 








Use the common information among 
different statuses of the dividing 
progress of targets. 


Accelerate the analysis of excitations 
from multiple angles. 


How do this? 


Analysis of scattering from 


Analysis of scattering from 
targets in free space. 


targets embedded in planar- 
layered media. 


GPU-Acceleration 
(Under Development) 








Y Discretization of targets 
dividing into parts using the 
RWG basis functions; 

Y Discretization of the connection 
between components parts 
using mono-polar RWG basis 
functions; 

Y Combined-Field Integral 
Equation using Discontinuous 


Galerkin method for analyzing 
Multi-Scale Discretization; 

Y Multi-level CBFM for reducing 
the unknowns and accelerate the 
computation; 

Y Far-field scattering calculation. 


Y Discretization of targets dividing 
into parts using the RWG basis 
functions; 

Y Mixed Potential Integral Equation 
for analyzing the objects in planar- 
layered media using MoM; 

Y Discrete Complex Image Method 
for analyzing the Dyadic Green's 
Functions in layered media; 

Y Multi-level CBFM for reducing 
the unknowns and accelerate the 
computation; 

Y Far-field scattering calculation 
using reciprocal theorem. 


v Cuda, Culbas and 
OPENGL for acceleration 
using GPU. 








Analysis of scattering from 
targets in free space 


How do this? 
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Discretization of targets dividing into parts 
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Combined-Field Integral Eguation in Free Space 


Boundary 
Condition 


Scattering 
Field 


Integral 
Operators 
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Combined-Field Integral Eguation in Free Space 
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Calculation of Elements of EFIE-MoM 
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Calculation of Elements of CFIE-MoM 
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Singular 
Integral 





Integral-Eguation Discontinous Galerkin Method 
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Weak-Form 
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Integral Equation 
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Multi-level CBFM for Reducing Unknowns and Acceleration 


Multi-Level 


Grouping 





Original 


Expande 


First Level 











Second Level 


In CBFM, first we need to group the basis functions 
into different levels, then use the lower-level basis 
functions to generate the higher-level characteristic 
basis functions. 


To avoid the effect of truncation, an expansion as 
large as 0.1~0.2 wavelengths is used when generate 
the CBFs. 


Multi-level CBFM for Reducing Unknowns and Acceleration 


To gain the information of the lower-level basis 


functions in each group when illuminated by plane 

Planewave te i Na ih 
YT wave rom d erent angles, we need to excite the 
Y» basis functions using plane waves from a series of 


angles. 
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Multi-level CBFM for Reducing Unknowns and Acceleration 
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Analysis of scattering from targets in free space: Numerical Result 1 





The target is comprised of 5 parts, a box as its main body and four pyramids which are truncated at the middle of them. The 
frequency is 525 MHz. Specitically, the body of the target and the upper and lower parts of the pyramid are discretized with 
elements whose lengths are 0.125, 0.0875 and 0.055 of the wavelength, respectively. The five stages of the target are analyzed, 
where the pyramids are connected to the body in the first stage, break away trom the body one by one in the next four 
stages. The spatial basis functions are subdivided into 12 blocks for level-1 and 48 blocks for level-2 via the octree algorithm 
to construct the CBFs. As the target undergoes a transition from the second through fifth stages, the CBFs and the reduced 
matrices associated with the blocks that have experienced movements are updated, while the other CBFs and the reduced 
matrices are saved for reusing. 


Analysis of scattering from targets in free space: Numerical Result 1 
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G8 — Polarized Mono-static RCS of Example 1 Obtained By Using CBFM #g— Polarized Mono-static RCS of Example | Obtained By Using CBFM 


Analysis of scattering from targets in free space: Numerical Result 1 


EE From the results, we can find 

| [| Mehod | i | xz | w |: that both the results calculated 
MON 237 10. by using 2L-CBFM and 1L-CBFM 

agree well with the reference for 
both polarizations. The CPU 
time, number of unknowns and 


the memory usage are reduced 
significantly. 





Analysis of scattering from targets in free space: Numerical Result 2 





For the second example, a PEC propeller in free space is analyzed. The length, width and thickness of the two blades are 10m, 
0.4m and 0.1m, respectively. The length and radius of the shaft are 1.6m and 0.15m, respectively. A box is used as the base, the 
length, width and height of which are 2.5m, 1.5m and 0.3m, respectively. The trequency is 480 MHz. The propeller is truncated 
into two parts at z = 1m, and their three stages are analyzed. The propeller is shown in Figure 5. The three stages correspond 
to the blades rotating at the angles 0°, 22.5? and 45°, respectively, while the base remains stationary. The spatial basis 
functions are subdivided into 58 blocks for level-1, and 2/ blocks for level-2 via the octree algorithm utilized in the CBFM. For 
the second and third stages, the CBFs and the reduced matrices are associated with the blocks that undergo changes because 
of the rotation are updated, while the others are saved so that they can be reused. 


Analysis of scattering from targets in free space: Numerical Result 2 
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Analysis of scattering from targets in free space: Numerical Result 2 


6690 3334 682.9 169.6 


MOM 


682.9 
247: 9347.7 





From the results, we can find 
that both the results calculated 
by using 2L-CBFM and 1L-CBFM 
agree well with the reference for 
both polarizations. The CPU 
time, number of unknowns and 
the memory usage are reduced 
significantly. 


Summary 


[1] Peterson AF, Ray SL, Mittra R. Computational methods for electromagnetics. New York (NY): IEEE Press; 1998. 

[2] Prakash VVS, Mittra R. Characteristic basis function method: a new technique for efficient solution of method of moments matrix equations. Microwave Opt Technol Lett. 
2003;36(2):95-100. 

[3] Lucente E, Monorchio A, Mittra R. An iteration-free MoM approach based on excitation independent characteristic basis functions for solving large multiscale electromagnetic 
scattering problems. IEEE Trans Antennas Propag. 2008;56(4):999-100/. 

[4] Mittra R. Computational electromagnetics. New York (NY); Springer; 2003, p. 209. 

[5] Laviada J, Las-Heras F, Pino MR, et al. Solution of electrically large problems with multilevel characteristic basis functions. IEEE Trans Antennas Propag. 2009;57(10):3189-3198. 

[6] Delgado C, Manuel Felipe C, Mittra R. Efficient multilevel approach for the generation of characteristic basis functions for large scatters. IEEE Trans Antennas Propag. 2008;56(7):2134— 
ET. 

[/] Laviada J, Gutierrez-Meana J, Pino MR, et al. Analysis of partial modifications on electrically large bodies via characteristic basis functions. IEEE Antennas Wirel Propag Lett. 
2010;9(1):834—837. 

[8] Fenni I, Roussel H, Darces M, et al. Application of the characteristic basis function method (CBFM) on a non-uniform mesh to the solution of large-size dielectric scattering problems. 
European Conference on Antennas and Propagation. Hague: IEEE; 2014; p. 2423-2426. 

[9] Nguyen H, Roussel H, Tabbara W. A coherent model of forest scattering and SAR imaging in the VHF and UHF-band. IEEE Trans Geosci Remote Sensing. 2006;44(4):838—848. 980 Y. 
SU ET AL. 

[10] Peng Z, Wang XC, Lee JF. Integral equation based domain decomposition method for solving electromagnetic wave scattering from non-penetrable objects. IEEE Trans Antennas 
PrOoaGg.2011,59(9) 3326-3336. 

[11] Peng Z, Lim KH, Lee JF. A discontinuous Galerkin surface integral equation method for electromagnetic wave scattering from nonpenetrable targets. IEEE Trans Antennas Propag. 
2013;61(7):3617-3628. 

[12] Cockburn B, Karniadakis GE, Shu CW. Discontinuous Galerkin methods: theory, computation and applications. Berlin: Springer; 2000. 

[13] Arnold DN. An interior penalty finite element method with discontinuous elements. SIAM J Numer Anal. 1982;19(4):742-760. 

[14] Houston P, Perugia |. Mixed discontinuous Galerkin approximation of the Maxwell operator. SIAM J Numer Anal. 2004;42(1):434—459. 

[15] Zhao Y, Ding D, Chen R. A discontinuous Galerkin time-domain integral equation method for electromagnetic scattering from PEC objects. IEEE Trans Antennas Propag. 
2016;6A(6):2410—2417. 

[16] Ubeda E, Rius JM. Novel monopolar MFIE MoM-discretization for the scattering analysis of small objects. IEEE Trans Antennas Propag. 2006;54(1):50—57. 

[17] Vechinski DA, Rao SM. A stable procedure to calculate the transient scattering by conducting surfaces of arbitrary shape. IEEE Trans Antennas Propag. 1992;40(6):661-665. 

[18] Zhao K, Vouvakis MN, Lee JF. The adaptive cross approximation algorithm for accelerated method of moments computations of EMC problems. IEEE Trans Electromagn Compat. 
2005;47(4):763-713. 

[19] Maaskant R, Mittra R, Tijhuis A. Fast analysis of large antenna arrays using the characteristic basis function method and the adaptive cross approximation algorithm. IEEE Trans 
Antennas Propag. 2008;56(11):3440—3451. 


